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High-performance, integrated spectrometers operating in the far-infrared and submillimeter ranges 
promise to be powerful tools for the exploration of the epochs of reionization and initial galaxy formation. 
These devices, using high-efficiency superconducting transmission lines, can achieve the performance of 
a meter-scale grating spectrometer in an instrument implemented on a 4 inch silicon wafer. Such a de- 
vice, when combined with a cryogenic telescope in space, provides an enabling capability for studies of the 
early universe. Here, the optical design process for Micro-Spec (p-Spec) is presented, with particular at- 
tention given to its two-dimensional diffractive region, where the light of different wavelengths is focused 
on the different detectors. The method is based on the stigmatization and minimization of the light path 
function in this bounded region, which results in an optimized geometrical configuration. A point design 
with an efficiency of ~90% has been developed for initial demonstration and can serve as the basis for 
future instruments. Design variations on this implementation are also discussed, which can lead to lower 
efficiencies due to diffractive losses in the multimode region. © 2014 Optical Society of America 
OCIS codes: (050.1940) Diffraction; (120.6200) Spectrometers and spectroscopic instrumentation; 

(220.2740) Geometric optical design; (300.6340) Spectroscopy, infrared. 
http://dx.doi.org/10.1364/AO.53.001094 


1. Introduction 

Far-infrared (IR) and submillimeter (15 pm-1 mm) 
spectroscopy provides a powerful tool to probe a wide 
range of environments in the universe. In the past 30 
years, many space-based observatories, such as Spit- 
zer, Herschel, and the Infrared Space Observatory 
(ISO) [1—3], have opened the far-IR window to the 
universe, revealing rich line and continuum spectra 
from objects ranging from interplanetary dust par- 
ticles to major galactic mergers and young galaxies 
in the early universe. Discoveries made by these 
observatories have provided unique insights into 
physical processes leading to the evolution of the 
universe and its contents. However, many questions 
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remain unanswered regarding the very first mo- 
ments of the universe, as well as galactic, stellar, 
and planetary formation. 

Micro-Spec (p-Spec) is proposed as a technology 
concept to enable new discoveries in the far-IR 
spectral range. p-Spec will be a high-sensitivity, 
direct-detection spectrometer operating in the 450- 
1000 pm regime. The light from the telescope is 
coupled to the spectrometer chip via a broadband slot 
antenna (Fig. 1) and directed to the power splitters 
and delay network realized in microstrip transmis- 
sion lines. These microstrip elements are connected 
to a planar waveguide multimode region via transi- 
tions, which also serve as internal feed apertures. 
The planar waveguide multimode region has two 
internal antenna arrays, one for transmitting and 
one for receiving the spectrum as a function 
of wavelength. Absorber structures lining the 
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Fig. 1. Layout of the |i-Spec module. The power is coupled into the instrument using a broadband antenna (left). It is then transmitted 
through a low -loss superconducting transmission line to a divider and a phase delay network, which creates a retardation across the input 
to the multimode region (in light blue). The feed horns will radiate a converging circular wave, which will concentrate the power along the 
focal surface, with different wavelengths at different locations. The outputs are connected to a bank of order-sorting filters to disentangle 
the various orders. 


multimode region are used to control the power emit- 
ted into large angles or reflected from the receiver 
antenna array. Finally, each receive antenna in the 
multimode region is coupled to a bank of order- 
sorting filters terminated in microwave kinetic 
inductance detectors (MKIDs) for readout. The in- 
strument, with a resolution 1Z ~ 1200, will be inte- 
grated on a 4 inch diameter (100 mm) silicon chip. 
This reduction in size is largely realized through 
the use of single-mode microstrip delay lines, which 
can be relatively compactly routed on the silicon wa- 
fer. Relative to counterparts realized in free space, 
the required transmission line lengths are more com- 
pact by a factor of the medium’s effective index. For 
all these reasons, p-Spec can become an important 
instrument capability under the low background con- 
ditions realized in balloon and space-borne platforms 
such as the Space Infrared Telescope for Cosmology 
and Astrophysics (SPICA) [4], The spectral range 
of the implementation studied here is set by the 
materials employed for the transmission line and 
detectors. Application of this approach is presently 
confined to wavelengths 1 > 250 pm by the gap fre- 
quency of available low-loss superconductors. These 
include niobium (Nb) and niobium-titanium nitride 
(NbTiN) for the transmission line structures and 
molybdenum nitride (MoN) for the detectors. 

Two distinct electromagnetic modes are used in the 
p-Spec transmission spectrometer: microstrip and 
parallel-plate waveguide [5] (see Fig. 4, right). Micro- 
strip is a single-mode transmission line structure that 
is utilized in its dominant quasi-transverse electro- 
magnetic mode (TEM) and in this context is analogous 
to single-mode fiber optic. It is used to realize the elec- 
trical delay line network. The parallel-plate wave- 
guide has a TEM symmetry, and the region serves 
as a multimode spatial power combiner. Its function 
can be compared to a grating spectrometer [6], in 
which a plane wave is reflected from the grating and 
the phase of each partial wave scattered from the 
rulings is a linear function of position across the gra- 
ting. p-Spec thus differs from similar technologies by 


the order of processing of the light in the spectrometer. 
For instance, in a Rowland spectrometer the required 
phase retardation is generated by reflection from 
the grating grooves [6—10] , whereas in Z-Spec, which 
is an example of planar Rowland grating architecture, 
propagation occurs in parallel-plate waveguides 
[11—14]. A last comparison can be made with one- 
dimensional bootlace lenses found in microwave 
practice [15-18], which p-Spec builds upon for sub- 
millimeter wave applications. 

In this work, the detailed optical design of the 
p-Spec transmission spectrometer is presented. 
We begin in Section 2 with an outline of the multi- 
mode power-combiner region optical design, then de- 
scribe the internal two-dimensional feed structures 
used to couple from microstrip lines to this multi- 
mode region in Section 3, and finish by bringing 
these elements together to compute the modeled re- 
sponse of the entire antenna array in Section 4. 

2. Detailed Spectrometer Design 

The performance of the spectrometer is optimized by 
the selection of the x and y coordinates of each radiator 
and a differential phase shift on each delay line. In the 
case of the Rowland spectrometer, two points on the 
focal surface can be chosen and wavelengths assigned 
to them [11,12] . The x andy positions of the grooves in 
the grating can be selected to provide zero phase er- 
rors at each of these stigmatic points. A similar con- 
structive approach can be used in the case of p-Spec, 
where each emitter has an additional degree of free- 
dom, the phase shift at the radiator, which allows gen- 
erating a third stigmatic point. In the case of the 
Rowland grating, the selection of wavelengths at 
the two stigmatic points determines the dispersion 
of the grating. 

In p-Spec, the selection of the wavelengths of three 
stigmatic points determines both the dispersion 
and the linearity of the dispersion of the system. A 
constraint on the selection is required if linear 
dispersion is desired. In the case presented here, 
the stigmatic points and their associated wavelength 
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were chosen to provide relatively linear dispersion 
and good imaging over the wavelength range of 
interest. The positions and wavelengths of the stig- 
matic points are free parameters and can be selected 
to optimize aspects of the instrument performance. 
A reasonable choice, which is adopted here, is to se- 
lect points that minimize the overall root-mean- 
square (RMS) phase error for all spectrometer 
outputs (Fig. 3). 

In the prototype version of the instrument de- 
scribed here, in order to fit four spectrometers in a 
100 mm silicon wafer, each in a 25 mm x 25 mm re- 
gion, the maximum radius allowed is R — 1.25 cm. 
This choice is not optimal but enables several designs 
to be experimentally studied within a single mask set 
and the development of the required fabrication proc- 
ess. The overall design was therefore carried out with 
the parameters shown in Table 1, where values for an 
initial design (A) and its optimized version (B) here 
adopted as the baseline are shown. 

With these parameters, the design procedure de- 
scribed was implemented in MATLAB through the 
goal attainment method of Gembicki [19]. This 
method was chosen because its implementation ena- 
bles the resolution of vectorial equations subject to 
nonlinear constraints that arise in the case presented 
in this paper. The following set of equations describes 
the propagation of light through the spectrometer, 

F m (x) = F m^ rn = 1 3 N e , (1) 

where F m represents the electrical delay that must be 
equal to F^, an integer multiple of the wavelengths, at 
the three stigmatic points. N e represents the number 
of radiators, and 3 N e equations are needed to calcu- 
late the electrical delay of all radiators at each of 
the three stigmatic points. The vector of the un- 
knowns, x, consists of the coordinates X; and y, of 
the N e emitters’ centers and the electrical path 
lengths in silicon, R®, for each feed horn delay line. 


Table 1. Spectrometer Parameter Summary* 


Global Design Parameters: 
Minimum frequency, fi = c/X 0 \ 


450 

[GHz] 

Maximum frequency, f 2 = c/X o2 


650 

[GHz] 

Silicon, relative permittivity, e r 


11.7 


Maximum wavelength in silicon, 


195 

[pm] 

A = -hi / a/ e rl l r 

Minimum wavelength in silicon, 


135 

[pm] 

h = 

Blaze point wavelength in silicon, 


164 

[pm] 

Xb = a/ X\X% 
Resolving power, 1Z 


65 


Grating Order, M 


1 


Design Variation: 

A 

B 


Number of emitters, N e 

65 

65 


Number of receivers, N r 

47 

65 


Antenna array element spacing 

179 

161 

[pm] 

or pitch, p 

Multimode region, physical radius, R 

1.25 

1.05 

[cm] 


“Speed of light in vacuum, c = 3.0 x 10 8 m/s. 


Constraints are necessary to limit the search of the 
solution in the correct parameter space. In fact, each 
of the system’s equations is of order 2 and the search of 
the solution was restricted within X \/ 40 from the geo- 
metric optics solution for the grating solution (the 2 R 
dotted arc on the left side of Fig. 2). Using a numerical 
as opposed to an analytical solution thus allows us to 
analyze and incorporate the channel isolation and 
higher-order diffractive response into the optimiza- 
tion (see Section 3). Indeed, while the feed phase cen- 
ters reside on the position of the analytical solution, it 
will be shown that the array response is not uniform 
in angle. 

All quantities in Eq. (1) are defined as follows: 


F{x) = 


Rf+Rli-Rl 
Re i +R e 2i- Re 2 

Rl+Rm-Rs 


i = 1 N e 


is the light path function matrix (3 N e x 1), whose first 
2 N e rows are associated with stigmatic points 1 and 
2 respectively, and the last N e rows with the blaze 
point B. All quantities in the F(x) matrix represent 
electrical lengths in the transmission medium. For 
example, the terms in the first N e rows associated 
with stigmatic point 1 are indicated below: 

• R\ is the electrical path length of the i-th silicon 
microstrip delay line; 

• R\i = [V (R sin /A -Xi ) 2 + (R cos /h -yd 2 } x yU/G 
is the electrical path length of the i-th feed horn 
from stigmatic point 1 in the silicon parallel-plate 
waveguide multimode region; 


y 



left side three radiators can be seen, which point to the blaze point, 
B. The path from each radiator’s phase center to the first and sec- 
ond stigmatic points are also indicated by solid lines for the array’s 
i-th and center reference feed. 
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• R\ — [R y/( 1 + sin f) i) 2 + cos 2 /fi] x ^fFpT r is the 
reference electrical path length of stigmatic point 1 
from the central feed horn (i — (N e + 1) /2), in the sil- 
icon multimode region. 

A 0.45 pm thick monocrystalline silicon dielectric 
substrate is employed on the chip. The electrical 
length in the parallel-plate waveguide region is com- 
puted from the product of the index of refraction in the 
medium and the physical length defined by the geom- 
etry in Fig. 2. For example, R e u — A 1; ^/e,.p r , where e r 
and fi r « 1 are the effective relative permittivity and 
permeability specifying the propagation in the guid- 
ing medium. The parallel-plate waveguide conductor 
and ground plane are well approximated and modeled 
as perfect electric conductors (PEC). A surface react- 
ance was assigned to the superconductor (niobium) to 
account for the kinetic inductance and model its influ- 
ence on the impedance and propagation constant in 
the microstrip transmission line structures. The 
superconducting material was modeled as a 0.2 pm 
thick layer with a London penetration depth X L — 
95 nm [20]. Similar definitions apply to stigmatic 
point 2 and the blaze point B. Finally, 


F* = 


( i_ 

0 - 

( f - 


¥ )ma o1 
^)m o2 
^ )ma oB 


i = 1, ...,N e , 


where M — 1 represents the order of the grating 
chosen for this design and A 0 the free-space wave- 
lengths of the three stigmatic points. Here, the pres- 
ence of subscript o signifies the radiation wavelength 
in free space and without indicates in the dielectric 
medium as defined in Table 1. In order to maximize 
the instrument efficiency, the emitters are tilted by 
an angle a so as to point to the center of the focal sur- 
face (Fig. 2): 


R e i + \ y /(R-x i ) 2 +y?-R £ 


X £ r lA r 


i- Ne ^ V \MA oB . 


( 4 ) 


Now, the terms y (R - xj 2 + yf express the distance of 
the i-th feed horn from the blaze point. Because the 
emitters lie on the 2 R circle, this quantity is exactly 
equal to 2 R. Also R B is equal to 2 R given the axisym- 
metric configuration of the instrument and the 
particular location chosen for the blaze point. There- 
fore, these two terms cancel out and one readily 
obtains 


Re = MA oB [i-^^\, (5) 

which is a linear relation between the electric length 
in silicon, R®, and the feed horn number, i. 

To find the precise spacings of the emitters along 
the grating circle needed in order to satisfy the stig- 
matic equation for the pair of symmetric off-axis 
stigmatic points at the specified high and low wave- 
lengths, our numerical solution method shows that 
the spacing between adjacent antennas is not a con- 
stant but grows quadratically as a function of posi- 
tion. In design example A, a change of ~13 pm 
from the nominal design value, p — 179 pm, is ob- 
served. This is the same behavior found by Rotman 
and Turner [15], who derived an analytical solution 
for the analogous design of a bootlace lens and found 
that the problem reduces to solving an ordinary 
quadratic equation, from which one can determine 
the antenna positions. 

As mentioned in the beginning of this section, the 
diffractive solution can be found by restricting the al- 
gorithm to searching for radiator positions within a 
limited distance from the ideal Rowland circle. The 
constraint reads 


= - arctan 


( y ‘ V 

V I x i I + Rj 


( 2 ) 


y /(x i -R) 2 +y 2 -{2R + z) 


< /fi/40. 


( 6 ) 


A particular solution is found when placing the 
emitter antennas on the 2 R circle, using an axisym- 
metric configuration, and placing the blaze point on 
the optical axis. This provides a partial diagonaliza- 
tion of the problem, as in [15]. Under these hypothe- 
ses, indeed, the last N e equations of Eq. (1), which 
refer to the blaze point, can be explicitly rewritten 
as shown below: 

\j (R sin ~ *i ) 2 + (R cos () B -yj 2 - R B 
x - N& 2 (3) 

where fi B is the angle between the y axis and the line 
from the center of the Rowland circle to the blaze point 
(Fig. 2). Here !> B — n/2, which yields 



where z is the phase center correction provided by 
finite -element simulations (z — 15 pm, see Section 4) 
and Ai/40 ~ 5 pm is chosen to provide the upper and 
lower bounds. This is equivalent to a phase error 
of (2n/A 0 i) x (4 o i/40) rad ~ 1/6 rad. 

Having selected the nine parameters, i.e., the coor- 
dinates of the stigmatic points and their wave- 
lengths, we proceed to determining the locations 
and phase shifts of the radiators. The geometry of 
a Rowland spectrograph was chosen as a starting 
point, with the radiators on a surface of radius 2 R 
and the detectors on a circle of radius R at a distance 
of 2 R from the center of the radiators. Finally, the 
spacing of the radiators is a free parameter and will 
typically be chosen to limit diffraction losses. We be- 
gin with a first radiator near the center of the radi- 
ator circle. Next, we choose an adjacent radiator and 
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perturb its position and differential phase, which 
results in arriving at each of the three stigmatic 
points. This procedure is continued for the radiators 
on each of the delay lines. This results in zero phase 
error at the stigmatic points by construction, low 
RMS phase errors between the stigmatic points, 
and diffraction-limited performance in the multi- 
mode region over the design spectral range. The 
RMS phase error, </? RMS , is illustrated in Fig. 3 for 
design A, and it was calculated by assuming a uni- 
form weight for each radiator at each point on the 
focal plane: 


<PRMs(0/) — 


-N„ 

E 


((Pij - {(p(Oj ))) 2 

N„ 


1/2 


(7) 


Here, ip^ is the relative phase of each transmitter, 
{<p(Oj)) — 0 is the relative phase of the central trans- 
mitter (this is zero by construction as the central ra- 
diator is used as a reference, i.e., F* N hlj/ , 2 = 0), N e is 
the number of transmitters, and 6j represents the an- 
gle corresponding to each of the points in which the 
focal plane was discretized. The phase between the 
i-th transmitter and the j-th point on the receiver 
focal plane is defined as cpij — 2n s /~e^jT r / X 0 j x F, and 
X 0 j is the ffee-space wavelength. 


3. Antenna Feed Response 

Conceptually, each microstrip feed structure is a two- 
dimensional analog of an adiabatic feed horn. To 
achieve the desired performance, the feed needs to 
match impedance, modal shape, and angular resolu- 
tion over the design frequency band. An adiabatic 
impedance taper [ 21 ] for the input microstrip line 
width to the output is achieved via a Hecken taper 
[22], which is a near-optimum matching section in 
the sense that it approaches an “equal-ripple” re- 
sponse without the use of abrupt discontinuities. 
These properties lower its sensitivity to fabrication 
variations and tolerances. This feed section serves 



Fig. 3. RMS phase error between the average and the phase of all 
other transmit feed paths evaluated as a function of receiver posi- 
tion with the design parameters in Table 1 for design A. The RMS 
phase error vanishes at the three stigmatic points by construction 
and remains low in between. This ensures diffraction-limited per- 
formance over the power combiner’s design spectral range. The 
feed positions for the receive array in this design configuration 
are indicated by circular symbols in the figure. Addition of the ra- 
diators indicated by crosses can be used to improve the throughput 
in future implementations. 


as a transition between microstrip lines and paral- 
lel-plate TEM waveguide regions. The latter is com- 
monly referred to as the parallel-plate guide, and the 
microstrip planar waveguide model [5] can be used to 
compute its propagation and impedance properties. 
The transformer ratio, effective dielectric constant, 
and desired return loss set the necessary length of 
the impedance taper and were synthesized from 
readily available analytical forms and the response 
numerically simulated. A diagram of the feed aper- 
ture and spacing geometry is provided in the insert 
of Fig. 4 (left). In illuminating the multimode dielec- 
tric-filled parallel-plate waveguide region, the micro- 
strip outputs play the role of the antenna apertures 
and need to be sized to minimize diffraction at the 
upper band edge. 

When used as a two-dimensional antenna element, 
the aperture geometry not only sets the radiation 
pattern but also plays a role in determining the 
return loss [23], This can be seen by the following 
consideration: The angular acceptance of the struc- 
ture is constant over the band since the wavelength 
in the medium is less than the effective aperture 
width, w. Considering the scaling of the reflection am- 
plitude as a function of feed size [ 24 ] , the desire for a 
low return loss at the low end of the band is in tension 
with the needs to limit diffraction ( w < a 0 / ^Je r p r ) 
and efficiently radiate (w > 4 0 /2 ^/e r p r ). This behavior 
effectively determines the smallest aperture size that 
can effectively radiate [25]. The reflection from the 
end of the microstrip taper scales as a power of the 
wavelength over the microstrip width in two dimen- 
sions and is larger than that arising from the imped- 
ance taper, which can readily be made subdominant 
by increasing the transformer’s length. Thus, to ad- 
dress the reflection from the aperture, the output 
geometry of the aperture was apodized by adding a 
constant radius flare [26]. The resulting radiator 
element return loss, >20 dB, is achieved over a 2 : 1 
bandwidth. Over the wavelength range where 
2 P <: 4, only the first-order diffraction occurs [27]. 
The influence of higher-order diffraction is mitigated 
through the feed apodization, the coupling geometry 
of arrays, and the use of absorbing boundaries on the 
multimode region walls. 


4. Antenna Array Response 

An efficient and accurate expression is desired for the 
evaluation of the structure’s response during numeri- 
cal optimization. To evaluate the antenna array’s re- 
sponse, a classical phased array analysis approach is 
used. The emitters and receivers in the multimode 
region can be considered as a pair of opposing two- 
dimensional antenna arrays. By considering the 
Poynting vector as a function of spatial dimension, 
d, the following expression for the electric field is 
adopted for the computation of the array’s response: 


N, 


E(r,0) = £ 


E 0 (6) 


^1 + (M)T -' 4 


. e /(kT-¥'y)_ 


( 8 ) 


1098 


APPLIED OPTICS / Vol. 53, No. 6 / 20 February 2014 



Parallel-Plate Waveguide 



0.2pm Nb 

Vacuum 

f{\ 





A 

,\ 

1 0.45pm Si 

1 1 


0.2pm Nb 


Microstrip Line 


Fig. 4. (Left) The computed feed horn angular response at 430 GHz. The response of an emitting feed (aperture size 73 pm) is evaluated in 
the far field (||r|| ~ 61 = 438 pm) and normalized to the magnitude of the E field at zero degrees. The feed array geometry is provided in the 
figure insert. The feed’s phase center is indicated by a black filled circle. (Right ) Cross section of the parallel-plate waveguide (upper) and 
microstrip (lower) transmission line geometries. The electric field is depicted with arrows to indicate the dominant modal symmetry. An 
adiabatic transition serves as a mode converter between these modes and defines the feed’s angular response in the multimode power 
combined region. 


Here, r is the separation distance vector between the 
center point of the i-th emitter and a point in the mul- 
timode region, E 0 is the E-field amplitude from the 
beam pattern angular distribution ( 9 ) at the center 
of the i-th feed horn (Fig. 4), w is the effective width 
of the feed horns, k = {2k Je r ^ r / A n ) ■ n, is the wave- 
number vector in the medium, and n, is the unity vec- 
tor normal to the i-th emitter. The phase cp^ includes 
the contributions arising from the feed aperture’s 
phase, the phase delay, and the reference to the cen- 
tral emitter. 

The field amplitude’s dependence on radial separa- 
tion is physically motivated by the behavior of the 
fields in transitioning from the near to the far field 
[28]. Close in to the feed aperture the fields are es- 
sentially constant, and in going to large distance 
the amplitude decreases with distance. Physically 
useful limiting cases to consider are cL — 3, which 
yields a scaling of the E-field inversely with distance 
for a spherical wave, and d — 1, which yields a scal- 
ing independent of distance for a lossless transmis- 
sion line. For d — 2, the case of primary interest 
here, it is important to note that the E-field ampli- 
tude scales as the square root of the separation dis- 
tance, ||r||, and is anticipated for a cylindrical 
geometry from the asymptotic form of the Hankel 
functions. This expression conserves power flow, is 
analytic over the parameter range of interest, and 
is a computationally efficient representation of the 
structure’s behavior during synthesis of the array re- 
sponse. The functional form and scaling indicated in 
Eq. (8) were validated by high-frequency structure 
simulation (HFSS) finite-element simulations of 
the structure’s response with distance. These simu- 
lations were used to derive the feed aperture’s 


effective width, w, and to determine the location of 
the phase center for radiation emitted by the 
structure. 

The power is then computed from 

P(r,9) = ±E*(r,0)-E(r,0), ®) 

where i] 0 — 3770 is the wave impedance of free space. 
Figure 5 illustrates the power distribution at 
450 GHz, as was calculated from Eqs. (9) and (8). 
On the left side, it can be seen that the power peaks 
at the first stigmatic point as required, and then it 
quickly drops to zero along the rest of the focal plane. 
At separations from the feed elements greater than 
~2 mm, the field transitions to the far field and the 
electric field amplitude is observed to asymptotically 
fall as ||r|| -1 / 2 . Simulations show the power is con- 
served throughout the multimode region. In design 
example A, the ratio of the power emitted by the feed 
horns to the power received by the antennas is equal 
to ~30%. The feed response that is not intercepted by 
the antenna arrays is terminated by the absorber 
structures on the sides of the multimode region. This 
also enables any residual reflected power to be termi- 
nated in a controlled manner. In addition, a second- 
order diffraction peak, visible in the lower lefthand 
quadrant near the emitters, is responsible for part 
of this loss. 

Improved performance can be achieved with differ- 
ent trade-offs between the antennas and the array 
geometry, which increase the throughput and reduce 
the effects of diffraction (see Table 1, design B). If the 
two stigmatic points are moved away from the blaze 
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x (cm) x (cm) 

Fig. 5. Power distribution in the multimode region at a frequency of 450 GHz. On the left (Design A), the first-order peak is visible at an 
angle 8 ss 10°, whereas the second-order diffraction peak shows up at 8 a -66° (8 is defined in Fig. 2). In the optimized version (design B) on 
the right, the higher-order diffraction peak has now disappeared and the throughput is as high as 87%. 


point on the optical axis along the Rowland circle 
by 4.2%, the instrument reaches its highest efficiency 
of 87% when the following applies: The radius is 
decreased by 16% (R — 1.05 cm), the number of 
receivers is increased by 18 (N r — 65), and the 
emitter pitch, p, is set to 161 pm. The resulting in- 
crease in the receiver array acceptance angle is bet- 
ter matched to capture the transmitted beams given 
the apodization function of the radiators. This is the 


primary reason for the multimode coupling-efficiency 
increase observed in this example. As an aside, this 
change in array geometry reduces the magnitude of 
the second-order diffraction peak and moves it out of 
the range of angles of interest. For additional detail, 
see Fig. 5, rightmost inset. 

Figure 6(a) shows the emitter antenna array’s nor- 
malized power response evaluated at the three 
stigmatic wavelengths for the optimized multimode 




1 10 20 30 40 50 60 


Receiver number 


(b) 


Fig. 6. Design B angular response, (a) The emitter antenna array’s angular power response is normalized by the total radiated power and 
plotted for the three stigmatic wavelengths. The peak directivity Dq(A,8) is the angular response of the array relative to an isotropic 
radiator and is indicated on the plot, (b) Power isolation between the receiver channels computed for X D ~ 560 pm, where the power in 
each channel is normalized to the power in the central channel. 
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Table 2. 

Computed Coupling Efficiency (2 p > 

A) 

Signal Frequency 

Guide Wavelength 

Coupling Efficiency 
Design A Design B 

[GHz] 

[pm] 

[-] 

H 

450 

X x = 195 

0.30 

0.91 

541 

X B = 164 

0.30 

0.87 

650 

i 2 = 135 

0.34 

0.55 


region design. Table 2 summarizes the coupling 
efficiency between receiver and transmitter antenna 
arrays for the two design variations considered in 
this section. The interchannel power isolation is pre- 
sented in Fig. 6(b) for the center channel of the opti- 
mized design. The observed finite isolation, <0.48, 
arises from the design goal to fully sample the 
spectrum and maximize the structure’s optical 
throughput. The resulting signal correlation can be 
addressed in calibration of the spectrometer re- 
sponse. Finally, the RMS phase error now reaches 
a maximum of 0.011 at/7 — 75.8° and 104.3° and van- 
ishes at the three stigmatic points (/! = 67.8°, 90° 
and 112.2°). 

5. Conclusions 

A theory and design methodology has been described 
for p-Spec, a superconducting microstrip-based 
direct-detection spectrometer, in which the phase 
retardation is generated by propagation through 
transmission lines of different length. A proof-of- 
concept design was realized to demonstrate the fea- 
sibility of the instrument. The location of the three 
stigmatic points was chosen to minimize the overall 
RMS phase error for all spectrometer outputs, 
whereas the location of both radiators and receivers 
was found by setting the light path function equal 
to integer multiples of the three wavelengths associ- 
ated with the three stigmatic points. The structure’s 
response was evaluated through an E-field model 
validated by finite-element simulations. The results 
show that, through a geometric optimization of the in- 
strument, a maximum throughput of ~90% can be 
achieved in the multimode region. The attributes of 
this design approach, namely reduced spectrometer 
envelope and high throughput, have the potential 
to enable a new class of space-born instrumentation 
that will reveal new details of the early universe. 
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